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Propagation of balance-sheet or cash-flow insolvency across financial institutions may be 
modeled as a cascade process on a network representing their mutual exposures. We derive 
rigorous asymptotic results for the magnitude of contagion in a large financial network and 
give an analytical expression for the asymptotic fraction of defaults, in terms of network 
characteristics. Our results extend previous studies on contagion in random graphs to 
inhomogeneous directed graphs with a given degree sequence and arbitrary distribution 
of weights. We introduce a criterion for the resilience of a large financial network to the 
insolvency of a small group of financial institutions and quantify how contagion amplifies 
small shocks to the network. Our results emphasize the role played by "contagious links" 
and show that institutions which contribute most to network instability in case of default 
have both large connectivity and a large fraction of contagious links. The asymptotic results 
show good agreement with simulations for networks with realistic sizes. 
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1 Introduction 



The recent financial crisis has highlighted the complex nature of linkages between financial 
institutions. From balance-sheet exposures, to more opaque obligations related to over-the- 
counter derivatives like credit default swaps, such linkages propagated and amplified financial 
distress. Initial losses in one asset class -mortgage backed securities- turned into losses that 
threatened the stability of the whole financial system. More than 370 of the almost 8000 US 
bank companies have failed since 2007. This was clearly an episode of large default contagion, 
if we compare to a number of 30 defaults in the period 2000 — 2004, and no defaults occurred 
in the period 2005 - 2006. 

The acknowledgement of different types of interbank connections and the associated conta- 
gion mechanisms led to an increased advocacy to account for network effects when discussing 
regulatory requirements |24l 1X4] . In the growing body of work dedicated to systemic risk, sev- 
eral distress propagation mechanisms have been pinpointed, including primarily two types of 
insolvency: balance-sheet insolvency and cash-flow insolvency. A bank is said to be balace- 
sheet insolvent if the value of its liabilities exceeds the value of its assets and, it is said to be 
cash-flow insolvent if it cannot meet its contractual payment obligations arrived at maturity. 

Cascades of insolvencies can be understood as domino effects and are a type of potent 
financial contagion that is quantifiable. 

Most investigated, balance-sheet insolvency contagion can be described as follows: party A 
has a balance sheet exposure to party B, where we understand by exposure the maximum loss 
incurred by A on its balance sheet claims upon the default of B. If B defaults, the capital of 
party A must absorb the corresponding loss. If the capital cannot withstand the loss, party A 
becomes balance sheet insolvent . 

Another contagion mechanism that can be described in similar terms is represented by 
cascades of cash-flow insolvencies. Over-the-counter derivatives markets are prone to such 
type of cascades. Indeed, in these markets parties deal directly with one another rather than 
passing through an exchange. As such, they are subject to the risk that the other party 
does not fulfill its payment obligations. Consider two parties A and B, such that A has a 
receivable from party B upon the realization of some event. If B does not dispose of enough 
liquid reserves, it will default on the payment. Now consider that B has entered an off-setting 
contract with another party C, hedging its exposure to the random event. If C is cash-flow 
solvent, then the payment will flow through the intermediary B and reach A. However, if C 
is cash-flow insolvent and defaults, then the intermediary B might become cash-flow insolvent 
if it depends on receivables from C to meet its payment obligations to A. As it turns out, 
the length of such chains of intermediaries in certain over-the-counter markets, like the credit 
default swap market, is significant [151 H6], thereby increasing the probability of a cascade of 
cash-flow insolvencies 

Some of previous work, mostly in the economics and sociology literature, investigates cas- 
cades on networks in a generic context; relevant references include Morris [42] , Kleinberg |34| . 
Jackson and Yariv [25] and Watts |50j . These models consider, in one form or another, a 
mechanism by which an agent decides to adopt one of two states depending on the state of its 

1 The trigger of a cash-flow insolvency cascade in a chain of intermediaries may be a counterparty that is not 
contractually required to post collateral when all other entities have this obligation. 
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neighbors and a threshold which measures its susceptibility to this direct influence. Propaga- 
tion of insolvency in banking networks fall under the irreversible version of this model |50| [33] 
; default is not reversible, unlike the case of agents playing a network game who can revise 
their decisions [32] . The default threshold of a given bank depends on its level of capital (resp. 
liquidity reserve), the state of balance sheet (resp. cash-flow) solvency of its direct counterpar- 
ties and the linkages to them. More recent work studies network formation as a result of the 
interplay between benefits from creating links and negative network externalities [35]. 

In the finance literature, in the context of banking systems, contagion effects and network 
externalities have been investigated in both in theoretical [U [TUl \19\ [22], 135] and empirical 
studies |171 120] 139] 133] . Network externalities are -implicitly or explicitly- present in various 
early discussions of systemic risk (see e.g., Hellwig [25J, Kiyotaki and Moore [33J, Rochet and 
Tirole [36]) through the interlinkages between balance sheets. Allen and Gale [1] pioneered 
the use of network models in the study of the stability of a system of interconnected financial 
institutions. Their results were extended in various directions by Lagunoff and Schreft [35] and 
Leitner [37]. 

One of the central problems tackled in this literature is understanding the relation between 
the cascading behavior of the network and the underlying topology: Is the network such that 
the state of a small number of nodes will propagate to a large number of nodes, or will contagion 
die out quickly? 

The discussion was either dominated by highly stylized networks, whose structure turns out 
to be quite different from the heterogenous structure of real networks -financial networks are 
particularly heterogenous, as many of the empirical studies [131 H3 E] make a clear case- or 
the results were heuristic in nature and based on mean field approximations |50l 122"] . Notable 
exceptions come from the random graphs literature, where cascade models are investigated on 
graphs with given degree sequences [2 [38] • One crucial aspect that does not appear in the 
previous literature is the heterogeneity of weights: whether interpreted as exposures or receiv- 
ables, these linkages carry weights with a heavy tailed distribution. This point -corroborated 
by simulations |17j- prevents from reducing the analysis of contagion in banking networks to 
the case where a node's aggregate exposure is distributed equally across counterparties as in 

[221132]. 

In light of insights coming from empirical studies and simulations, we redefine the problem 
as relating the cascading behavior of financial networks both to the local properties of the nodes 
and to the underlying topology of the network. Since balance-sheet and cash-flow insolvency 
cascades are similar from a mathematical modeling point of view, the challenge lies not so 
much in analyzing a model that is flexible enough to represent both these types of insolvency 
cascades, but in proposing a model that can mimic the empirical properties of these different 
types of networks and that is tractable enough to be able to prove theorems about the cascading 
behavior. 

Our problem is set form the point of view of a regulator who observes the network. Our 
primary goal is not to identify nodes posing the highest systemic risk. Clearly, when one 
knows the entire network and assuming the network is fairly small, those could be identified 
by extensive simulations. Alternatively, for threshold models of contagion in large networks 
one could use approximation algorithms to find sets of most influential nodes [32, 43j. Indeed, 
problems like the Influence Maximization Problem have been shown NP-hard to approximate 
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within a factor 1 — 1/e + e for all e > [32] , but under certain sub- modularity conditions there 
exists a greedy 1 — 1/e — e - approximation algorithm for this problem |43j . 

The fundamental question that we tackle, unanswered so for realistic networks, is how to 
identify the features that make nodes systemically important. The obvious purpose is the need 
to set rules that would mitigate such features and consequently systemic risk. 

Our approach is to consider an ensemble of networks in which one can prescribe each 
node's connectivity and characteristics that are relevant to the respective cascade mechanisms: 
balance sheet insolvency cascades depend on capital ratios and the asset side of the balance 
sheets; cash-flow cascades depend on liquidity reserves and cash flows related to positions in 
the trading book. When the number of banks is large, cascades on networks belonging to this 
ensemble behave in a way dictated by the prescribed characteristics. 

1.1 Summary 

In this paper we develop techniques for analyzing default cascades in random weighted directed 
networks with arbitrary degree sequences, in which a set of local features can be prescribed for 
each node in the network. 

Our contribution is to derive rigorous asymptotic results for the magnitude of contagion in 
such networks and give an analytical expression for the asymptotic fraction of defaults. Our 
results apply to a wide variety of topologies and provide analytical insights into the nature of the 
relation between network structure, local characteristics of nodes and contagion in large-scale 
networks. 

For simplicity, we formulate our results in terms of balace-sheet insolvency cascades in a 
network of financial institutions with interlinked balance sheets, where losses flow into the 
asset side of the balance sheets. Similar techniques may be used for analyzing cascades of 
cash-flow insolvency in over-the-counter markets, as briefly discussed above and detailed in 
|40| . From now on, we refer to balance-sheet insolvency simply as insolvency. 

Our proof is based on a coupling argument: We construct a related multigraph -a weighted 
configuration model- which leads to the same number of defaults as in the original contagion 
process but is easier to study because of its independence properties. The contagion process in 
this model may then be described by a Markov chain. Generalizing the differential equation 
method of Wormald [51 J to the case where the dimension of the Markov chain depends on the 
size of the network we show that, as the network size increases, the rescaled Markov chain 
converges in probability to a limit described by a system of ordinary differential equations, 
which can be solved in closed form. This enables us to obtain analytical results on the final 
fraction of defaults in the network. 

These results generalize previous ones on diffusions in random graphs with prescribed degree 
sequence to the case of inhomogeneous and weighted random directed graphs with arbitrary 
degree sequences. Related problems are the problem of existence of a giant component in 
random graphs [THIE] and bootstrap percolation problem. Bootstrap percolation process^] is a 

2 A bootstrap percolation process on a graph G is an "infection" process which evolves in rounds. Initially, 
there is a subset of infected nodes and in each subsequent round each uninfected node which has at least r 
infected neighbors becomes infected and remains so forever (The parameter r > 2 is fixed.). 
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very simple models of difusions which have been studied on a variety of graphs, such as trees [8], 
grids [261 E] > hypercubes [6] , as well as on several distributions of random graphs [31 EH [9] . 

Another important result of our work -and probably the most important from the regulatory 
point of view- is the introduction of a measure of resilience of a financial network to small initial 
shocks. The contribution of each node to systemic risk is quantifiable in terms of its connectivity 
and local characteristics. Our measure may be used as a tool for stress testing the resilience of 
interbank networks in a decentralized way [1J and as an assessment tool of the capital adequacy 
of each bank with respect to its exposures. 



1.2 Outline 

The paper is structured as follows. Section [2] introduces a model for a network of financial 
institutions and describes a mechanism for default contagion in such a network. Section [3j gives 



our main result on the asymptotic magnitude of contagion. Section 3.4 uses this result to define 
a measure of resilience for a financial network: We show that when this indicator of resilience 
crosses a threshold, small initial shocks to the network -in the form of the exogenous default 
of a small set of nodes- may generate a large-scale cascade of failures, a signature of systemic 
risk. Section [4] illustrates, through concrete examples, how the resilience measure allows us to 
quantify and predict the outcome of contagion on one sample network generated from a random 
network model that mimics the properties of a real interbank exposure network analyzed in 
|17j . We observe that networks with the same average connectivity may amplify initial shocks 
in very different manners and their resilience to contagion can vastly differ. In particular, the 
relation between 'connectivity' and 'contagion' is not monotonous. Technical proofs are given 
in Appendix [A} 



2 A network model of default contagion 

In this section, we first introduce a model of a financial network, then describe the default 
cascade on this network, and finally the probabilistic setting we use throughout the paper. 



2.1 Counterparty networks 

Interlinkages across balance sheets of financial institutions may be modeled by a weighted 
directed graph G = (V,e) on the vertex set V = {1, . . . , n} = [n], whose elements represent 
financial institutions. The exposure matrix is given by e 6 ]R nxn , where the ij-th entry e(i,j) 
represents the exposure (in monetary units) of institution i to institution j. Table [l] displays 
a stylized balance sheet of a financial institution. The interbank assets of an institution i are 
given by 

A(i) :=J2<hj)- 
j 

Note that ^2je(j,i) represents the interbank liabilities of i. In addition to these interbank 
assets and liabilities, a bank may hold other assets and liabilities (such as deposits). 
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Assets 


Liabilities 


Interbank assets 


Interbank liabilities 


E,-e(i,j) 


Ej e U> *) 




Deposits 






Other 


Net worth 


assets 




x(i) 


c(i) = 7(i)A(t) 



Table 1: Stylized balance sheet of a bank. 



The net worth of the bank, given by its capital c(i), represents its capacity for absorbing 
losses while remaining solvent. We will refer to the ratio 



as the "capital ratio" of institution i, although technically it is the ratio of capital to interbank 
assets and not total assets. 

An institution is insolvent if its net worth is negative or zero, in which case we set 7(2) = 0. 

Definition 2.1 (Financial network). A financial network (e, 7) on the vertex set V = [n] is 
defined by 

• a matrix of exposures {e(i, j)}i<i,j<n, 

• a set of capital ratios {7(i)}i<j< ra - 

In this network, the in- degree of a node i is given by 

d-(i):=#{jeV\ e(j,i)>0}, 
which represents the number of nodes exposed to i, while its out-degree 

d + (i):=#{j£V\ e(i,j)>0} 

represents the number of institutions i is exposed to. 

The set of initially insolvent institutions is represented by 

D (e,7) = {^e^l 7(0 = 0}- 
The next section defines the default cascade triggered by nodes in Do( e )7)- 
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2.2 Default contagion 

In a network (e, 7) of counterparties, the default of one or several nodes may lead to the 
insolvency of other nodes, generating a cascade of defaults. 

Starting from the set of initially insolvent institutions Do( e > 7) which represent fundamental 
defaults, we define a contagion process as follows. 

Denoting by R(j) the recovery rate on the assets of j at default, the default of j induces 
a loss equal to (1 — R(j))e(i, j) for its counterparty i. If this loss exceeds the capital of i, 
then i becomes in turn insolvent. Recall that c(i) = j(i)A(i). The set of nodes which become 
insolvent due to their exposures to initial defaults is 

Di(e j7 ) = {*€F| -y(i)A(i)< £ (1 - R(j))e(i,j)}. 

This procedure may be iterated to define the default cascade initiated by a set of initial defaults. 

Definition 2.2 (Default cascade). Consider a financial network (e, 7) on the vertex set V = [n]. 
Set Do(e,7) = {i 6 V | = 0} of initially insolvent institutions. The increasing sequence 

(Bfc(e, 7), k > 1) of subsets of V defined by 

o k (e, 1 ) = {iev\ i(i)A(i)< Yl (i-^'M.;)} 

ieD fc _i(e,7) 

is called the default cascade initiated by Do( e >7)- 

Thus Dfc(e,7) represents the set of institutions whose capital is insufficient to absorb losses 
due to defaults of institutions in Dfc_i(e, 7). 

It is easy to see that, in a network of size n, the cascade ends after at most n — 1 iterations. 
Hence, D n _i(e, 7) represents the set of all nodes which become insolvent starting from the 
initial set of defaults H>o(e, 7). 

Definition 2.3. Consider a financial network (e,7) on the vertex set V = [n]. The fraction of 
defaults in the network (e,7) (initiated by Do( e >7)) is given by 

n 

The recovery rates R{i) may be exogenous or, as in Eisenberg and Noe |19j . determined 
endogenously by redistributing assets of a defaulted entity among debtors, proportionally to 
their outstanding debt. As noted in }48tll7j. the latter scenario is too optimistic since in practice 
liquidation takes time and assets may depreciate in value due to fire sales during liquidation. 
As argued in [171 120], when examining the short term consequences of default, the most realistic 
assumption on recovery rates is zero: Assets held with a defaulted counterparty are frozen until 
liquidation takes place, a process which can in practice take months to terminate. 

For simplicity, we assume from now on that recovery rates are constant for all institutions: 

R(i) = R, Vt G V. 
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2.3 A random network model 

Empirical studies on interbank exposures |13} [T7] show such networks to have a complex and 
heterogeneous structure characterized by heavy-tailed (cross-sectional) distributions of degrees 
and exposures. 

Given a description of the large-scale structure of the network in statistical terms, it is 
natural to model the network as a random graph whose statistical properties correspond to 
these observations. 

Consider a sequence (e n ,7„) n >i of financial networks, indexed by the number of nodes n, 
where d+ = {d^(i)}f =i (resp. d~ = {d~(i)}f =1 ) represents the sequence of in-degrees (resp. 
out-degrees) of nodes in e n . We now construct a random network E n such that e n may be 
considered as a typical sample of E n . 

Definition 2.4 (Random network ensemble). Let G n ( e n) be the set of all weighted directed 
graphs with degree sequence d+, d~ such that, for any node i, the set of exposures is given by 
the non-zero elements of line i in the exposure matrix e n . Let (S7, A, P) be a probability space. 
We define E n : Q — > Q n(On) <H random directed graph uniformly distributed on Q n {e n ). 

We endow the nodes in E n with the capital ratios ~f n . Then for all i = 1, . . . , n, 
{E n (i,j), E n (i,j)^0} = {e n {i,j), e n (i,j)^0} P - a.s. 

#{j G V, E n (j, i)>0} = d+(j), and #{j G V, E n (i,j) ? 0} = d~(i). 

Definition |2.4| is equivalent to the representation of the financial system by an unweighted 
graph chosen uniformly among all graphs with the degree sequence (d^, d~), in which we assign 
to the links emanating from node i the set of weights {e n (i,j) > 0}. 

3 Asymptotic results 

We consider a sequence of random financial networks as introduced above. Our goal is to 
study the behavior of a n (E n ,7„) which represents the final fraction of defaults in the cascade 
generated by the set of initially insolvent institutions, i.e., Do(E n ,7 n ) = {i G [n] \ ^y n (i) = 0}. 

3.1 Some probability-theoretic notation 

We let No be the set of non-negative integers, i.e., No = N U {0}. For non-negative sequences 
x n and y n , we describe their relative order of magnitude using Landau's o(.) and 0(.) notation. 
We write x n = 0(y n ) if there exist N G N and C > such that x n < Cy n for all n > N, and 
x n = o(y n ), if x n /y n -)■ 0, as n -)■ oo. 

Let {X n } ng N be a sequence of real-valued random variables on a sequence of probability 
spaces {(r2„, P n )} ng N. If c G K is a constant, we write X n A c to denote that X n converges in 
probability to c. That is, for any e > 0, we have P n (|X n — c\ > e) — > as n — > oo. 

Let {a n } n£ N be a sequence of real numbers that tends to infinity as n — > oo. We write 
X n = o p (a n ), if \X n \/a n converges to in probability. Additionally, we write X n = O p (a n ), 
to denote that for any positive- valued function u>(n) — > oo, as n — > oo, we have PdA^/a^ > 
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uj(n)) = o(l). If £ n is a measurable subset of Q n , for any n G N, we say that the sequence 
{£ n } ng N occurs wi/i /ii<//i probability (w.h.p.) if P(£ n ) = 1 — o(l), as n — )• 00. 

Also we denote by Bin(/c,p) denotes a binomially distributed random variable correspond- 
ing to the number of successes of a sequence of k independent Bernoulli trials each having 
probability of success equal to p. 

3.2 Assumptions 

Consider a sequence (e n ,7 n ) n >i of financial networks, indexed by the number of nodes n. Let 
m n denote the total number of links in the network e n : 

n n 

m n := ^2d+(i) = 

i=i i=i 

The empirical distribution of the degrees is defined by 

MnC/jfc) : = e [n] I =j,d~(i) = k}. 

n 

From now on, we assume that the degree sequences d+ and d~ satisfy the following condi- 
tions analogous to the ones introduced in |41j . 

Assumption 3.1. For each n G N, d+ = {{d£(i))f =1 } and d~ = {(^n(*))I=il are sequences 
of non- negative integers with XrfLi^nW = SiLi^n(O) an d such that for some probability 
distribution \i on Nq independent of n and with finite mean A := ^ ■ fc A;) = ■ fc fc) G 
(0,oo), the following holds: 

1. fjinU-i ^) ~~ * f° r every j, > as n — > 00 ; 

2- E[Li(4W) 2 + K«) 2 = 0(n). 

Note that, in particular the second assumption implies (by uniform integr ability) that 
m n /n — > A, as n — > 00. 

We now present our assumptions on the exposures. Let us denote by S n (i) the set of 
all permutations of the counterparties of i in the network e ra , i.e., permutations of the set 
{j G [n] I e n (i,j) > 0}. For the purpose of studying contagion, the role of exposures and 
capital ratios may be expressed in terms of default thresholds for each node. 

Definition 3.2 (Default threshold). For a node i and permutation r n G S n (i) which specifies 
the order in which z's counterparties default, the default threshold 

n k 

Qn{hT n ) ■■= min{/c > I j n (i) e n (i, j) < }Xl - R)e n (i, T n (j))} (1) 

i=i i=i 

measures how many counterparty defaults i can tolerate before it becomes insolvent (in the 
financial network (e n ,j n )), if its counterparties default in the order specified by r n . 
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We also define 

(■ U a\ _ #{(h T n) I i ^ H Tn ^ ^{i), (t IkM = h = fc ' Q n(i,T n ) = 9} , . 



We will see in Section A. 2 that for n large, p n (j,k,9) gives the fraction of nodes with degree 
(j, k) which have a default threshold equal to 9, in the random financial network E n . 
In particular for 9 = 1, 

n^ n (j,k)jp n (j,k,l) 

is the number of exposures of nodes with degree (j, k) which exceed the capital of the exposed 
node, i.e., exposures which in case of default of the initial node always lead to the insolvency 



of the exposed node. These links play a crucial role (as we will see in Section 3.4) and we call 
them contagious links. 

Definition 3.3 (Contagious link). We call a link contagious if it represents an exposure 

larger than the capital of the exposed node: 

n 

(1 - R)e n (i,j) > c n (i) = 7n{i)^2e n (i,j). 

j'=i 

From now on, we assume that p n (j, k, 9) has a limit when n — > oo. 
Assumption 3.4. There exists a function p : Nq — > [0, 1] such that for all j, k, 9 £ No (9 < j) 

p n (j, k, 9) -> p(j, k, 9), as n -> oo. 



Some examples of exposures for which this assumption is fulfilled are given in Section [4. 1| 
Under this assumption, we will see in Appendix [A] that p(j, k, 9) is also the limit in probability 
of the fraction of nodes with degree (j, k) which become insolvent after 9 of their counterparties 
default. In particular, 

• p{j, k, 0) represents the proportion of initially insolvent nodes with degree (j,k); 

• p(ji 1) represents the proportion of nodes with degree (j, k) which are 'vulnerable', i.e., 
may become insolvent due to the default of a single counterparty. 

3.3 The asymptotic magnitude of contagion 



Consider a sequence (e n ,7„) n >i of financial networks satisfying Assumptions 3.1 and 3.4, and 
let (E n ,7 n ) n >i be their corresponding sequence of random financial networks, see Definition 
2.4[ Let us denote by 



PV, 7T, 9) := P(Bin(j, tt) > 9) = ]T Q tt'(1 - vr) 



i>e 

the distribution function of a binomial random variable Bin(j, tt) with parameters j and tt. 
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We define the function / : [0, 1] — > [0, 1] as 

J W := £ fc, vr, 0). ( 3 ) 

i,fc /X e=o 

Indeed, I(tt) has the following interpretation (when the network size goes to infinity): If the 
end node of a randomly chosen edge defaults with probability tt, I(tt) is the expected fraction 
of counterparty defaults after one iteration of the cascade. 
Let tt* be the smallest fixed point of I in [0, 1], i.e., 

tt* = inf{vr G [0,1] | I (it) = tt}. 

(The value tt* represents the probability that an edge taken at random ends in a defaulted 
node, at the end of the contagion process.) 

Remark 3.5. / admits at least one fixed point. Indeed, I is a continuous increasing function 
and, 

j,k 8=0 

since ^2Qp(j,k,9) < 1 by definition. Moreover, 

J(0)=^^^ p(i) fc,0)>0. 

j,k 

So the function / has at least a fixed point in [0, 1]. 

We can now announce our main theorem. 
Theorem 3.6. Consider a sequence (e n ,1n)n>i of financial networks satisfying Assumptions 



3.1 and 3.4 • and the corresponding sequence of random matrices (E n ) n >j defined on (O,^.,] 



as in Definition 2.4 Let tt* be the smallest fixed point of I in [0, 1]. 

1. If tt* = 1, i.e., if I(tt) > tt for all tt G [0, 1), then asymptotically almost all nodes default 
during the cascades 

a n (E n ,7„) 4 1. 

2. If tt* < 1 and furthermore tt* is a stable fixed point of I, i.e., I 1 (it*) < 1, then the 
asymptotic fraction of defaults is given by 

j 

a n (E n , 7n ) 4 »ti, k) vr*, 9). 

j,k 9=0 

A proof of this theorem is given in Appendix |A} 
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3.4 Resilience to contagion 



The resilience of a network to small shocks is a global property of the network which depends 
on its detailed structure. However, the above results allow us to introduce a rather simple and 
easy to compute indicator for the resilience of a network to small shocks. Consider a sequence 
( e ni7n)n>i of financial networks satisfying Assumptions 3.1 and 3.4 



Definition 3.7 (Network resilience). We define the network resilience function as 

I-EtM^W^,!) G (-oo,1]. 



The following result, which is a consequence of Theorem 3.6, shows that this indicator 
measures the resilience of a network to the initial default of a small fraction e of the nodes: 

Proposition 3.8. Consider a sequence (e n ,7 n ) ra >i of financial networks satisfying Assump- 
tions 3.1 and 3.4, and let (E n3 7 fl ) n >i be their corresponding sequence of random financial 



networks, see DefinitionWJX If 



i-^i(^Mj^,i)>o, 



j.k 



A 



(4) 



then for every e > there exists N e and p £ such that, if the initial fraction of defaults is smaller 
than p e then the final fraction of defaults is negligible with high probability: 



Vn > N F 



(a n (E n ,7 n ) <e)>l-£. 



Proof. Consider p be defined as the fraction of fundamental defaults 

P ■= ^2v(j,k)p{j,k,0). 

j,k 

We have 

I(a) = ££^£pO\M))90\M). 

j,k ' 8=0 

Using a first order expansion of j3(j, a, 9) in a at (when a — > 0), we obtain 

I3{j,a,9) = l{ d =o} + o?l{0=i} + o{a). 

Thus, 



1(a) = £ ^^(p( /■ k. 0) + ,,,>(./. /,-.!)) + of n ). 



3,k 



A 



Let a* be the smallest fixed point of 1(a). Given Condition Q, for a > and small enough, 
we have 

lim 1(a) = a (—II A p(j^ x) + o(a) < a, 



A 
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where we used the fact that if the fraction of fundamental defaults tends to zero, so does the 
fraction of out-going links belonging to fundamentally defaulted nodes. On the other hand we 
have seen that 1(0) > 0. Thus linip^o ol* = 0. 

Let us now fix e > 0. By continuity of the function g defined by 



9(a) ■■= ^2/J,(j,k)^2p(j,k,0)P(j,a,9), 

j,k 8=0 



appearing in Theorem 3.6, there exists p £ such that g(a*) < e/2 as soon as p < p £ . By Theorem 
3.6 we have that there exists an integer N £ such that, for n > N £ , 



P(K(E n , 7 „)- 5 (a*)| <e/2) > 1 - e, 
which completes the proof. □ 

The proof points out what is a natural and intuitive fact: If contagion does not spread to 
nodes with contagion threshold 1, then it will not spread at all. 

Theorem 3.9. Consider a sequence (e n ,1n)n>i of financial networks satisfying Assumptions 



3.1 and 3.4. and let (E„,7 n ) n >i be their corresponding sequence of random financial networks. 



i-£4^,o,m)<o, ( 5) 

j,k 

then with high probability there exists a set of nodes representing a positive fraction of the finan- 
cial system, strongly interlinked by contagious links (i.e., there is a directed path of contagious 
links from any node to another in the component), such that any node belonging to this set can 
trigger the default of all nodes in the set. 

Given the network topology, Condition Q sets limits on the fraction of contagious links 
p n (j,k,l), i.e., on the magnitude of exposures relative to capital. A proof of this theorem 
is given in Appendix |A| However, Condition [4] may be justified using the following heuristic 
argument. 

Remark 3.10 (Branching process approximation). We describe an approximation of the local 
structure of the graph by a branching process, the children being the in-coming neighbors: the 
root 4> with probability ^(k^) := ]Cj Mil ^) has an in-degree equal to k^. Each of these k^ 

vertices with probability has degree (j, k), and with probability equal to p(j, k, 1) default 

when their parent defaults. Let y be the extinction probability, given by the smallest solution 
of 

y = J2^- P (j,k,l)y k . (6) 

If Yljk \ pU> k, 1) < 1, then the smallest solution of (6) is y = 1 (the population dies out 
with probability one), whereas if 

\r- Mj, k)jk . 

2^ — i — PU. k, 1) > 1, 



j-k 



13 



there is a unique solution with y G (0, 1). 

Remark 3.11 (Too interconnected to fail?). We suppose that the resilience condition Q is 
satisfied. Let tt* be the smallest fixed point of I in [0, 1], when a fraction s of all nodes represent 
fundamental defaults, i.e., p(j, k, 0) = e for all j, k. 

We obtain then, by a first order approximation of the function /, that 

K = 7 u 1" 

By a first order approximation of the function 7r — > £V fc A:) ^2l =0 p(j, k, 9)(3(j, it, 6) 
giving the asymptotic fraction of defaults in Theorem |3.6[ we obtain that, for any p there 
exists e p and n p such that for all e < e p and n > n p 

P(|«n(E n , 7 n) - e(l + ^ fc) fc )| < p) > 1 - P- (7) 

Suppose now that initially insolvent fraction involves only nodes with degree (d + ,d~), and 
we denote 7r*(d + , d~) the smallest fixed point of I in [0, 1] in the case where p(d + , d~, 0) = e 
and p(j, fc, 0) = for all (j, k) ^ (d + , d~). Then we obtain that, for any p there exists e p and 
n p such that for all e < e p and n > n p , 

P |«n(E n , 7 n) - + — ^ fc ^i^i '/ - )[ < p > 1 - p. (8) 



A i_ Ej . fc « p(j)M) 

This simple expression shows that there are basically two factors that determine how small 
initial shocks are amplified by the financial network: the interconnectedness of the initial 
default -represented by its in-degree d~- and the susceptibility of the network, given by the 
factor Y,j,k a 1 )" 

4 Numerical results on finite networks 

The results of Section [3] hold in the limit of large network size. In order to assess whether 
these results still hold for networks whose size is large but finite, we now compare our theo- 
retical results with numerical simulations for networks with realistic sizes. In particular, we 
investigate the effect of heterogeneity in network structure and the relation between resilience 
and connectivity. We begin by presenting some examples of models for counterparty networks 



which satisfy Assumption 3.4 



4.1 Examples of networks which satisfy Assumption 3.4 



In this section we give two important examples of exposures which satisfy Assumption 3.4 



Example 4.1 (Independent exposures). Assume that for all n, the exposures of all nodes 
i G [n] with the same degree (j,k), i.e., 

{e n (i,l) > | = j,d~(i) = k) , 
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are independent and identically distributed (i.i.d.) random variables, with a law given by 
Fx(j,k), depending on j and k but independent of n. We assume the same for the sequence 
of capital ratios, i.e., {"f n (i) \ d^(i) = j, d~{i) = k} are i.i.d. random variables with a law given 
by Fj(j, k) which may depend on (j, k), but not on n. Then it is easy to see that, by the law 



of large numbers, Assumption 3.4 holds and the limit p(j, k, 9) is known (for all j, k, 0), 

j 0-1 

p(j, k, 0) = F(X e > 7 ~ I> - R ) X * ^ °)> 

i=i i=i 

where 7 is a random variable with law Fy(j,k), and (Xi)j =1 are i.i.d. random variables with 
law Fx(j, k) and independent of 7. 

Example 4.2 (Exchangeable exposures). Empirical observations of banking networks [TTIEJ 
|4"7] show that they are hierarchical, 'disassortative' networks [39J, with a few large and highly 
interconnected dealer banks and many small banks, connected predominantly to dealer banks. 
This can be modeled in a stylized way by partitioning the set of nodes into two sets, a set V 
of n D dealer banks, and a set Af of n N non-dealer banks. 

We assume that the exposures {e n (i, I) > | i G V}, and {e n (i,l) > | i G A"} are 
restrictions corresponding to the first (respectively ) elements of infinite sequences of 
exchangeable variables, where and denote the total number of exposures belonging 
to dealer and respectively non-dealer banks. Similarly, the capital ratios {j n (i) \ i G 2?} and 
{7 n (i) I i G A/"} are restrictions to the first n D (respectively n N ) elements of the sequence, 
independent of the sequence of exposures. 

We can extend this example to a finite number of classes of nodes represented by their 
degrees, and also drop the assumption of independence between exposures and capital ratios. 
We assume that within each class, the sequence of a node's exposures and capital ratios are 
exchangeable random variables. 

For each node i with d^(i) = j, d~(i) = k, we let 

Y n {%) := ({e n (i,j) > 0},7 n (z')) 

be a multivariate random variable with state space 9^*' C M. We assume that the law of 
the finite sequence 

{Y n (i) I ie [n],d+{i)=j,d-{i) = k} 

is invariant under permutation. 

Then the family {Y n (i) \ i G [n],d^(i) = j,d~(i) = k} 0< j k<M represents a family of finite 
multi-exchangeable systems, as defined by Graham |23) . It is proved in [23] that the conditional 
law of a finite multiclass system, given the value of the vector of the empirical measures of its 
classes, corresponds to independent uniform orderings of the samples within each class, and that 
a family of such systems converges in law if and only if the corresponding empirical measure 
vectors converge in law. 

Let us consider the empirical measure sequence 



A j,k . = 1 { d+(i)=j, d-(i)=k}$Y n (i) 
"/';.!./'•/••; 



0<j,k<M 
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We suppose that the family {Y n (i) | % G [n], <i^(i) = j, d„ (i) = k} 0< j k<M converges in law 
when n — > oo to an infinite multi-exchangeable system 

Km | i g [nU+(i) =J,d^(i) = 4o<i, fc <M = W' fc I Z > ^o^^m" (9) 

By |23l Theorem 2], the empirical measure converges in law to 

lim {A^ fc ) k (a^) . (10) 

n->oo I J 0<j,k<M I J 0<j,fc<M 

For an arbitrary Z G we define the function 

#{r\re^(j),e(Z,r)=9} 



h(Z,9) 



Thus, by Equation (10) giving the convergence of empirical measures and, the fact that the 



function h is bounded, we have 

Pn (j, k, 9) = E A ™ (h(Z, 6)) "-=r IE AJ '" (h(Z, 0)) = p(j, k, 9), 

with Z a random element of Q 5 - 7 ' and E n and E ' denoting expectation under the measures 
A^ k and A^'* respectively. A last observation is that Equation (bj) is verified in our two tiered 
example since the sequences of exposures in the network of size n are restrictions of infinite 
exchangeable sequences. 

4.2 Relevance of asymptotics 

Interbank networks in developed countries may contain several thousands of nodes. The Federal 
Deposit Insurance Corporation insured 7969 institutions as of 3/18/2010, while the European 
Central Bank reports 8350 monetary financial institutions in the Euro zone (80% credit in- 
stitutions and 20% money market funds). To assess the relevance of asymptotic formulae for 
studying contagion in networks with such sizes, we generate a scale- free network of 10000 nodes 
with Pareto distributed exposures using the random graph model introduced by Blanchard |llj . 
which can be seen as a static version of the preferential attachment model. In this model, given 
the sequence of out-degrees, an arbitrary out-going edge is assigned to an end-node % with 
probability proportional to the power d^(i) a where a > 0. This leads to positive correlation 
between in-degrees and out-degrees. 

The distribution of the out-degree in this model is a Pareto law with tail exponent 7 + : 

^(i) == #{* I d+(i) = j} /x+(j) ~ f + +\ 
and the conditional limit law of the in-degree is a Poisson distribution 

P{d -= k \d + =j) = e-^^f, 

with A(j) = j;((D+) t ^) ' wrier e D + denotes a random variable with law . The main theorem 
in [11] states that the marginal distribution of the in-degree has a Pareto tail with exponent 
7~ = 3—, provided 1 < a < 7 + . 
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In the sequel we investigate cascades on these finite networks. A crucial ingredient is the 
capital ratio of each bank. Throughout this section, the capital ratio is assumed to be bounded 
from below by a minimal capital ratio and, we consider the worst case scenario where all 
nodes have a capital ratio equal to the minimal capital ratio: 7CO = Tminj 

Vi 6 V. 

The distribution of this simulated network's degrees and exposures is given in Figure [l] and 
is based on the empirical analysis of the Brazilian network in June 2007 [T7] . 



In-Degree Dislribution 




Out-Degtee Distribution 




Exposure Size Distribution 




Figure 1: (a) The distribution of out-degree has a Pareto tail with exponent 2.19, (b) The 
distribution of the in-degree has a Pareto tail with exponent 1.98, (c) The distribution of the 
exposures (tail-exponent 2.61). 

On one hand we make a simulation of the default contagion starting with a random set 
of defaults representing 0.1% of all nodes (chosen uniformly among all nodes). On the other 
hand we plug the empirical distribution of the degrees and the fraction of contagious links 
into Equation Q for the amplification of a small number of initial defaults. Figure [2] plots 
the amplification for varying values of the minimal capital ratios. We find a good agreement 
between the theoretical and the simulated default amplification ratios. We can clearly see that 
for minimal capital ratios 7 m j n less than the critical value 7^ n , the amplification ratio increases 
dramatically. 

Figure [3] plots the simulated fraction of defaults in a scale free network, starting from the 
initial default of a single node, as a function of the in-degree of the defaulting node, versus the 
theoretical slope given in Equation (8). Recall that the theoretical slope Ylj k \ k p(ji ^> -Q 
in Equation (TsT) measures the susceptibility of the network to initial defaults. 



17 



Amplification of initial defaults 



25 - 



Heterogenous network - Simulation 

A Heterogenous network - Theoretical 



20 - 



0.3 

Y 



Figure 2: Amplification of the number of defaults in a Scale- Free Network. The in- and 
out-degree of the scale- free network are Pareto distributed with tail coefficients 2.19 and 1.98 
respectively, the exposures arc Pareto distributed with tail coefficient 2.6b n — 10000. 
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The agreement is particularly good for nodes with large in-degrees, pointing to large am- 
plification if such nodes were to default in a network with high susceptibility. 




Figure 3: Number of defaulted nodes 



4.3 The impact of heterogeneity 

In the previous examples we can compute the minimal capital ratio 7j^ ira above which the 
network is resilient under Condition Q. Two factors contribute to the sum in Condition Q: 
connectivity of the node, and its fraction of contagious links. We compare in Figure |4j the 
ratio by which contagion amplifies the number of initial default in three cases: a scale free 
network with heterogeneous weights (exposures) , a scale free network with equal weights and a 
'homogeneous' random network (the Erdds-Renyi random graph where every directed edge is 
present with a fixed probability) with equal weights. All three networks are parameterized to 
have the same average degree, i.e., the same total number of links. It is interesting to note that 
in our example the most heterogeneous network is also the least resilient, as opposed to the 
homogeneous Erdds-Renyi network with the same distribution of exposures. These simulations 
corroborate the role played both by the network topology and the heterogeneity of weights. 



4.4 Average connectivity and contagion 

One of the recurrent questions in the economics literature regards the impact of connectivity 
on resilience to contagion: Is increased global connectivity posing a threat to financial system 
or does it allow for more efficient risk sharing? While an influential paper by Allen and Gale 
[1] finds that resilience increases with connectivity, Battiston et al. [10] exhibit different model 
settings where this relation is non-monotonous. Our results show nonetheless that in a model 
with heterogenous exposures, the average connectivity is too simple summary statistic of the 
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Heterogeneity and Systemic Risk 



Scale Free Heterogenous Exposures - Simulation 
Scale Free Heterogenous Exposures - Theoretical 
Scale Free Equal Exposures - Simulation 
Scale Free Equal Exposures - Theoretical 
- Erdos Renyi Equal Exposures - Simulation 
Erdos Renyi Equal Exposures - Theoretical 




0.35 0.4 0.45 0.5 



Figure 4: Amplification of the number of defaults in a Scale-Free Network (in- and out-degree 
of the scale-free network are Pareto distributed with tail coefficients 2.19 and 1.98 respectively, 
the exposures are Pareto distributed with tail coefficient 2.61), the same network with equal 
weights and an Erdos Renyi Network with equal exposures n = 10000. 



topology to be explain resilience. To see this let us consider a simple example and use the 
asymptotic formula for the amplification. 

For simplicity, consider networks in which nodes' exposures are equal and 1/3 < jmin < 1/2 
such that p n (j, k, 8) = l/j=i 2}- We consider three cases of degree distributions. 

First, let /i n (l,3) = fi n (2,3) = )U n (4, 3) = /U n (5, 3) = 1/4. The average connectivity in a 
network with this degree distribution is 3 and the resilience measure is equal to 1 /4. 

Second, let p, n (l,2) = 2/3, /2 n (4, 2) = 1/3. The average connectivity is 2 and the resilience 
measure is equal to 1/3. 

Last, we take /2 n (4, 4) = 1 i.e., a regular graph with degree 4. Here the average connectivity 
is 4 and the resilience measure is 1. 

In all three cases a network constructed with the empirical degree distribution is resilient 
w.h.p.. Nonetheless, we clearly observe that the resilience measure does not depend on the 
average connectivity in a monotonous way. 

While in the case of [10] this non-monotonicity is obtained by introducing an ad-hoc mech- 
anism of 'financial accelerators' on top of the network contagion effects, in our case it stems 
from an intrinsic trade-off between risk-sharing and contagion which is inherent in the model. 

These examples show that the resilience of a network cannot be simply assessed by exam- 
ining an aggregate measure of connectivity such as the average degree or the number of links, 
as sometimes naively suggested in the literature, but requires a closer examination of features 
such as the distribution of degrees and the structure of the subgraph of contagious links. 
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5 Conclusions 



In this paper we have analyzed distress propagation in a financial system where banks strongly 
differ in local features, and links between financial institutions are heterogenous in nature. 

We obtained an asymptotic expression for the size of a default cascade in a large network 
with prescribed characteristics -degree sequence and local features of the nodes- extending pre- 
vious results for homogeneous undirected random graphs to heterogeneous, weighted directed 
networks. 

Our asymptotic results were corroborated with a simulation study of contagion on a network 
with large but realistic size: on a given network, the spread of distress can be predicted by 
our measure of resilience. Our sample network has the same empirical properties as a real 
interbank network, e.g., the Brazilian one. As we vary the capital ratio, the point where 
the resilience measure becomes negative, or otherwise said where the theoretical amplification 
explodes, closely predicts the point where a large cascade would ensue. This illustrates how one, 
given an interbank network, could use the resilience measure and the theoretical amplification 
as a supervisory tool. This point is further developed in our follow-up paper on stress testing 



The crucial question in the context of macro prudential regulation of banking systems is 
how to identify and mitigate those features that make nodes systemically important. We have 
identified institutions acting as potential hubs for default contagion as those highly connected 
and with a large fraction of contagious links. One natural way to mitigate the systemic impact 
of these nodes is to set minimal capital requirements with respect to contagious links. This 
point of view is different from the current capital requirements as defined by the Basel II 
accords. Currently, minimal capital depends on the risk weighted sum of exposures, where the 
risk weights are given by the counterparty default probability. However, the default probability 
is computed by internal models and may not take into account knock-on effects. Our results 
suggest that, for financial stability, minimal ratios of capital should be set with respect to 
contagious exposures. 

While the insolvency contagion investigated in this paper has been mostly associated with 
balance-sheet contagion, our results may be applied to players in the over-the-counter deriva- 
tives markets. Contagion is carried in these markets through intermediaries with a large fraction 
of critical receivables defined similarly to contagious links, i.e., receivables on which the inter- 
mediary depends to meet its own payment obligations [16]. We argue that financial stability 
would be significantly enhanced by setting lower bounds on liquid reserves with respect to 
critical receivables, for those nodes which are counterparties to a large number of contracts. 

Last, whereas in this paper we were concerned mostly with a network where the node's 
characteristics were observable, our results apply to the particular case where connectivities and 
weights are sequences of (exchangeable) random variables with arbitrary correlation structure. 
This corresponds to a setting where the modeler cannot observe the sequence of exposures, 
but rather has some belief over their distribution. Such a perspective would have for example 
a market participant who models the default probability of its counterparties and takes into 
consideration network effets. 
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A Appendix: Proofs 



In this section we present the proofs of Theorem 3.6 and 3.8 We begin by introducing a 
weighted configuration model -a multigraph related to the financial network- which has the 
same asymptotic behavior as the random financial network. We then show by a coupling 
argument that the default cluster in the weighted configuration model can be constructed 
sequentially. The contagion process in this sequential model may then be described by a 
Markov chain. We then generalize the differential equation method of Wormald \§T\ to the case 
where the dimension of the Markov chain depends on size of the network and we show that, 
as the network size increases, the rescaled Markov chain converges in probability to a limit 
described by a system of ordinary differential equations. We solve these equations and obtain 
an analytical result on the final fraction of defaults in the network. Finally, we show that, if 
the resilience measure is negative, the skeleton of contagious links percolates. 



Consider a sequence (e n ,7 n ) n >i of financial networks satisfying Assumptions 3.1 and 3.4 



A.l Link with the configuration model 

A standard method for studying random graphs with prescribed degree sequence is to consider 
(see e.g., O [29] ) a related random multigraph with the same degree sequence, known as the 
configuration model [12], then condition on this multigraph being simple. The configuration 
model in the case of random directed graphs has been studied by Cooper and Frieze [18]. 
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1 \ "'d n (l) 1 \ ■■ d n {i) 



1 : c n (l) <■ : c n (i) 



1 * 1 * .. <^ 

2 2 



■d n (n) 



n : c n (n) 



1 * ..<*+(„) 



{e(l,-)>0} {e(v)>0} 



{e(n,0 >0} 



Figure 5: Configuration model 



Proceeding analogously, we introduce a multigraph with the same degrees and exposures as 
the network defined above, but which is easier to study because of the independence properties 
of the variables involved. Conditioned on being a simple graph, it has the same law as the 
random financial network defined above. 

Definition A.l (Configuration model). Given a set of nodes [n] = {1, . . . ,n} and a degree 
sequence (d+, d~), we associate to each node i two sets: H+(i) representing its out-going half- 
edges and H~(i) representing its in-coming half-edges, with \H+{i)\ = d^(i) and \H~(i)\ = 
d~{i). Let = H+(i) and H~ = |J,- H^(i). A configuration is a matching of with H~ . 
To each configuration we assign a graph. When an out-going half-edge of node i is matched 
with an in-coming half-edge of node j, a directed edge from i to j appears in the graph. The 
configuration model is the random directed multigraph G* (e n ) which is uniformly distributed 
across all configurations (Figure [5]). 

It is easy to see that, conditional on being a simple graph, G* (e n ) is uniformly distributed 
on Q n {e n ). Thus, the law of G* n (e n ) conditional on being a simple graph is the same as the law 
of E n . 

In particular any property that holds with high probability (with probability tending to 
1 as n — > oo) for the random multigraph G* (e n ), holds with high probability on the random 
network E n provided 

liminfP(G£(e n ) is simple) > 0. (11) 



In particular, the condition J2?=i(dn (*)) 2 + (d n (i)) 2 = 0(n) implies (11), see [30J. 
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Remark A. 2. Janson p30] has studied, in the case of undirected graphs, the probability of the 
random multigraph to be simple. One can adapt the proof to the directed case and show that 



the condition ^r=i(^rt W) 2 + (d n {i)) 2 = 0(n) implies (11). Indeed, in the non-directed case 



Janson [30] proves that when m n := Y^i=i d n {i) —> oo, (d n (i) is the degree of node i) one has 
P(G*(n, (dn(*))i) ^ simple) = exp A;; - £)(Ay - log(l + A;,-)) J + o(l), 

where for 1 < i,j < n; X ti := V^M&EMW^3 _ The proof of these results is based on 

J ill Ti 

counting vertices with at least one loop and, pairs of vertices with at least two edges between 
them, disregarding the number of parallel loops or edges. The same argument applies to the 
directed case, and one can show that when m n := Ya=i (*) = Ya=i ~~ * °°> then 

P(G* (e n ) is simple) = exp -- ^ A« - ^(A^ - log(l + Ay)) + o(l), 



where for 1 < i,j < n; Aj 



y 7 dj(i)dn (i)d£ (j)dn (j) 



One can observe that a uniform matching of half-edges can be obtained sequentially: choose 
an in-coming half-edge according to any rule (random or deterministic), and then choose the 
corresponding out-going half-edge uniformly over the unmatched out-going half-edges. The 
configuration model is thus particularly appropriate for the study of contagion, as we will see 
in the proofs, since we can restrict the matching process to choosing only in-coming half-edges 
entering defaulted nodes. In doing so, one constructs directly the contagion cluster in the 
random graph given by the configuration model and endowed with the sequence of capital 
ratios. 

Due to this property, it is easier to study contagion on G* n (e n ) under conditions on the 



degree sequence for the assumption above (11) to hold, then translate all results holding with 



high probability to the initial network E n defined in Definition 2.4 



A. 2 Coupling 

We are given the set of nodes [n] and their sequence of degrees (d+, d~). For each node i, we 
fix an indexing of its out-going and in-coming half-edges, ranging in = {1, ... ,<i+(i)} 

and [d^(i)] respectively. Furthermore, all out-going half-edges are given a global label in the 
range [1, . . . , m n ], with m n the total number of out-going (in-coming) half-edges. Similarly, all 
in-coming half-edges are given a global label in the range [m n ]. 

Recall that T, n (i) denotes the set of all permutations of the counterparties of i in the network 
e n . The set of weights on the links exiting node i is given by 

W n (i):={e n (i,j)>0}. (12) 

For the sequence of edge weights and capital ratios, (e n ,7„), we generate the random graph 
G n (e n ,7 n ), by the following algorithm: 
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1. For each node i, choose a permutation 6 S n (i) uniformly at random among all per- 
mutations of node i 7 s out-going half edges. 

2. Color all in-coming and out-going half-edges in black. Define the set of initially defaulted 
nodes 

V := (J {<}. 

*>7n(*)=0 

Set for all nodes i £ [n]\V ; c(i) = j n (i) J2j^[ n ] e n(hj)- 

3. At step k > 1, if the set of in-coming black half-edges belonging to nodes in T>k-i is 
empty, denote Vf the set T>k-i- Otherwise: 

(a) Choose among all in-coming black half-edges of the nodes in T>k-\ the in-coming 
half-edge with the lowest global label and color it in red. 

(b) Choose a node i with probability proportional to its number of black out-going half- 
edges and set 7r n (A;) = i. Let i have I — 1 out-going half-edges colored in red. Choose 
its r^(Z)-th out-going half-edge and color it in red. Let its weight be w. If the node 
i £ Pfc-i and (1 — R)w is larger than i's remaining capital then T>^ = T>k-i UW- 
Otherwise, the capital of node i becomes c(i) — (1 — R)w. 

(c) Match node i's T^(7)-th out-going half-edge to the in-coming half-edge selected at 



step (3a) to form an edge. 



4. Choose a random uniform matching of the remaining out-going half-edges and match 
them to the remaining in-coming half-edges in increasing order and color them all in red. 

Lemma A. 3. The random graph G n (e n ,j n ) has the same distribution as G* (e n ). Furthermore 
the set T>f at the end of the above algorithm is the final set of defaulted nodes in the graph 
G ? n (e n ,7 n ) (endowed with capital ratios j n ). 

Proof. The second claim is trivial. Let us prove the first claim. For a set A, we denote by 
the set of permutations of A. Let cr+ and a~ be the random permutations in Sr mn j, representing 
the order in which the above algorithm selects the in-coming / out-going edges. At step k of 
the above construction, in-coming half-edge with global label cr~(k) is matched to out-going 
half-edge with global label cr+(/c) to form an edge. The permutation <t+ is determined by 
the set of permutations (r*)j = i v . ^ n and the sequence ir n of size m n , representing the (ordered) 
sequence of nodes selected at Step 3b (or Step [4] when the set of in-coming black half-edges 



belonging to nodes in is empty - assume we choose sequentially uniformly at random) of 

the algorithm (each node i appears in sequence 7r n exactly d+(i) times). 

It is easy to see that c+ is a uniform permutation among all permutations in S[ mw i, since 
(7^)1=1,...^ are uniformly distributed and at each step of the algorithm we choose a node with 
probability proportional to its black out-going half-edges. On the other hand, the value of 
o~~(k) depends in a deterministic manner on 

(e n ,7 n ,o-+(l),...,o-+(£;-l)). 
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The out-going half-edge with global label j is matched with the in-coming half-edge with 
global label (a~ ° (o"n)~ )U)- ^ n order to prove our claim it is enough to prove that the 
permutation (a~ o (c^) -1 ) is uniformly distributed among all permutations of m n . Indeed, for 
an arbitrary permutation £ belonging to the set Sr m i, we have that 

P KO') = £~Vn 0'))l <{k) = CH*n(k)) for all k < j) = - 1 

m n j -\- 1 

Indeed, conditional on the knowledge of (<7+(l), . . . , cr+(j — 1)), cr~(j) is deterministic. Also, by 
conditioning on Vfc < j, o+{k) = ^(er" (k)), then f^OnC?)) G T := [m n ]\{<r+(l), . . . , cr+(j- 
1)}, of cardinal m„ — j + 1. In the above algorithm, cr^(j) has uniform law over T ■ Then the 
probability to choose i~ x {u~ (j)) is mn ^ j+1 - 

By the law of iterated expectations, we obtain that 

n°n ° {<?ir l = = p(4 = r 1 o a-) = 



This and the fact that the last step of the algorithm is a conditionally uniform match conclude 
the proof. □ 

We can find the final set of defaulted nodes T>f of the above algorithm in the following 
manner: once the permutation is chosen, assign to each node its corresponding threshold 



n (i) = Q n (i,Tn) as in Definition 3.2, and forget everything about (e n ,j n 



Definition A. 4. Denote by G n (d+, d~, # n ) the random graph resulting from the above algo- 
rithm, in which we replace Step |3b| of the algorithm by the fact that node i defaults the first 
time it has 8 n (i) out-going half-edges colored in red, i.e., at step 

inf-{> > 1, such that 9 n (i) = #{1 < I < k, ir n (l) = i}}. 

Corollary A. 5. The random graph G n (d+, d~, 9 n ) has the same law as the unweighted skeleton 
of G n (e n ,7 n ). 

Let N n (j, k, 9) denote the number of nodes with degree (j, k) and threshold 9 after choosing 
uniformly the random permutations r n in the above construction. 

Lemma A. 6. We have (as n — > oo) 

N n (j,k,e) v 



tJ-(j,k)p(j,k,8). 



Proof. For any node i with with degree (j, k), the probability that its default threshold G n (i, r* 
be equal to 9 is 

<■ m I Qn{i,r) = 9} 

v n (i,6) := 

Then we have 



N n (j,k,9)= Be(u n {i,9)), 

i, dt(i)=j, d„(i)=k 
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where Be(-) denotes a Bernoulli variable. By Assumption 3.4 we have 



E[N n (j, k, 6)/n] = Hn(j, k)p n (j, k, 9) n -T fi(J, k)p(j, k, 9), 

i „ TA7/-7 n\ i i d+U)=j, d-(i)=k 0)(1 - v n (i, 9)) 

and Var [N n (j, k, 9) n = y ' J ' w 5 -)• 0. 

n z 

Now it is easy to conclude the proof by Chebysev's inequality. 

□ 

A. 3 A Markov chain description of contagion dynamics 

In the previous section, we have replaced the description based on default rounds by an equiv- 
alent one based on successive bilateral interactions. By interaction we mean matching an 
in-coming edge with an out-going edge. At each step of the algorithm described in last section, 
we have one interaction only between two nodes (banks), yielding at most one default. This 
allows for a simpler Markov chain which leads to the same set of final defaults. 

We describe now the contagion process on the unweighted graph G n (d+, d~, 9 n ) with thresh- 
olds (9 n (i) = @ n (hT~n))i<i<n in terms of the dynamics of a Markov chain. 

At each iteration we partition the nodes according to their state of solvency, degree, thresh- 
old and number of defaulted neighbors. Let us define Sn ' ,l {t), the number of solvent banks 
with degree (j, k), default threshold 6 and I defaulted debtors before time t. We introduce the 
additional variables of interest: 

• D 3 n k ' d (t): the number of defaulted banks at time t with degree (j, k) and default threshold 
0, 

• D n (t): the number of defaulted banks at time t, 

• D~(t): the number of black in-coming edges belonging to defaulted banks, 
for which it is easy to see that the following identities hold: 

D^ d {t) = fJ, n (j,k) Pn (j,k,9)- &n k ' e '\t), 

0<l<6 

D~(t) = ]T kD^ 9 (t)-t, 

j,k,o<e<j 

D n (t) = £ D^\t). 

j,k,0<G<j 

Because at each step we color in red one out-going edge and the number of black out-going 
edges at time is m n , the number of black out-going edges at time t will be m n — t . 
By construction, Y n (t) = (S^ k,e '\t)j represents a Markov chain. Let (T ni , 



' j,k,o<i<e<j 

be its natural filtration. We define the operator A as 

x A y = max(x, y). 
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The length of the default cascade is given by 

T n = inf{0 < t < m n , D-(t) = 0} A m n , (13) 

The total number of defaults is given by D n (T n ), which represents the cardinal of the final set 
of defaulted nodes. 

Let us now describe the transition probabilities of the Markov chain. For t < T n , there are 
three possibilities for the partner B of an in-coming edge of a defaulted node A at time t + 1: 

1. B is in default, the next state is Y n (i + 1) = Y n (i). 

2. B is solvent, has degree (J, k) and default threshold 9 and this is the (Z + l)-th deleted 

out-going edge and Z + 1 < 6. The probability of this event is — — ■ The changes 

for the next state will be 

S^ l (t + l) = S^ e ' l (t)-l, 
S j,k,e,i+i( t + !) = s% k > e > l+1 {t) + 1. 

3. B is solvent, has degree (J, k) and default threshold 9 and this is the 6-th. deleted out-going 
edge. Then with probability ^~ 9+l ^ Sn W we have 

S j,k,e,o-i( t + 1 j = S j,k,e,o-i(j.} _ i 



Let At be the difference operator: AtY := Y(t + 1) — Y(t). We obtain the following equations 
for the expectation of Y n (i + 1), conditional on J- n> t, by averaging over the possible transitions: 



E 
E 



A t S^ e '°\T n , t 



A, qj,k,8,l \ x* 



m,. 



{j-l + l)sk kAl -\t) (j-l)S^ l (t) 



k,e,i, 



m r , 



(14) 



The initial condition is 



5J.M,/( ) = Nn (j : k} 0)x(| = o)l(0 < < j)- 



Remark A. 7. We are interested in the value of D n (T n ), with T n defined in (13). In case 
T n < m n , the Markov chain can still be well defined for t £ [T n ,m n ) by the same transition 
probabilities. However, after T n it will no longer be related to the contagion process and the 
value D~ (t) , representing for t <T n the number of in-coming half-edges belonging to defaulted 
banks, becomes negative. We consider from now on that the above transition probabilities hold 
for t < m n . 

We will show in the next section that the trajectory of these variables for t < T n is close to 



the solution of the deterministic differential equations suggested by equations (14) with high 
probability. 
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A. 4 A law of large numbers for the contagion process 

Define the following set of differential equations denoted by (DE): 

„• „j,k,0,O( \ 

(^*.«.»)'(t) = - , 

(sl , wy(T) _ o-i+i)^-'w_b--o^(r) t (DE) , 

A — r A — r 

with initial conditions 

^M-'(O) = /x(i,%(i,fc^)l(Z = O)l(O<0<i). 
Lemma A. 8. The system of differential equations (DE) admits the unique solution 

y(r) := (V' W (r)), 
in the interval < t < A, with 

JAB^ r) := M (j,*)p0'.*,«)(f)(l-x) i " J (x) ,1 {0<^- (15) 



j,k,o<i<e<j 



s J 



Proof. We denote by DE^ the set of differential equations defined above, restricted to j Ak < K 
and by b(K) the dimension of the restricted system. Since the derivatives of the functions 
[s : > ,k ' ' 1 (t)) . Ak<K0<l<6< . depend only on r and the same functions, by a standard result in the 

theory of ordinary differential equations |27[ Ch.2, Thm 11], there is an unique solution of DE^ 
in any domain of the type (— e, A) x R 7 with R a bounded subdomain of IR^^- 1 and e > 0. The 
solution of (DE) is defined to be the set of functions solving all the finite systems (DE^)x>i- 
We solve now the system DE. Let u = u(r) = —ln{\ — r). Then u(0) = —ln(X), u is 
strictly monotone and so is the inverse function r = r(u). We write the system of differential 
equations (DE) with respect to u: 

( s i> k > e > l )'(u) = U-l + l)J*M-\u)-(j-l)si> k > e > l (u). 



Then we have 

d 



(gj.M.i+igCj-t-iXu-upjh = (j _ n^W( u ) e (i-»-i)(«~w(o)) 
du 



and by induction, we find 



r=0 

By going back to r, we have 



j*M(u) = e -0-0(«-«(o)) ( J _ (l - e -("-u(0))\ s ^Ar {u{Q) y 



l-r 



^.'(r) = (l-I)^X;^(0)({_:) 

r=0 ^ ' 



J-r\ f r l _ r 



32 



Then, by using the initial conditions, we find 



•s- 



Ml {r) = n(j, k) P (j, k, 9) [ J (i - y-\-yi {e>0} 



□ 

A key idea is to approximate, following Wormald |51| . the Markov chain by the solution of 
a system of differential equations in the large network limit [51,31]. We summarize here the 
main result of |51j . 

For a set of variables Y 1 , Y b and for U C define the stopping time Tjj = TjjiY 1 , Y b ) 

inf{t > l,{t/n;Y l {t)/n,...,Y b (t)/n) £ U}. 

Lemma A. 9 ( Theorem 5.1. in |51j). Let b > 2 be an integer and consider a sequence of real 
valued random variables (0^(i)}i<K&)t>o and its natural filtration T n ^- Assume that there is 
a constant Co > such that \YJi(t)\ < Con for all n, t > and 1 < I < b. For all I > 1 let 
fi : M. b+1 —¥ M. be functions and assume that for some bounded connected open set U C M. b+1 
containing the closure of 

{(0, 2i, ...,z b ) : 3 n such that P(V 1 < I < b, 1^(0) = z { n) / 0}, 

the following three conditions are verified: 

1. (Boundedness). For some function /3(n) > 1 we have for all t < Tjj 

max|F^ + l)-^(t)|</3(n). 

l<i<0 

2. (Trend). There exists Ai(ra) = o(l) such that for 1 <l <b and t < Ty 

\E[Yji(t + 1) - Y l n {t)\F n , t ] - Mt/n, Y*(t)/n, ^(t)/n)| < Ai(n). 

3. (Lipschitz). The functions {fi)\<i<b are Lip schitz- continuous on U. 
Then the following conclusions hold: 

(a) For (0, z%, £&) G U, the system of differential equations 

^7 = fl(s,zi,-,zi), l = l,...,b, 

has a unique solution in U, z\ : K — > R, which passes through zi(0) = Z\, for I = 1, . . . , b, 
and which extends to points arbitrarily close to the boundary ofU. 

(b) Let A > Ai(n) with A = o(l). For a sufficiently large constant C, with probability 1 — 

0(^exp(-^))^e W 

sup (Yl{t) - nz l n {t/n)) = O(An), 

0<t<o-(n)n 
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where z n (t) = (z^(t), . . . ,z^(t)) is the solution of 

dz n 



dt 



f(t,z n (t)) z n (0) = Y n (0)/n 



and a(n) = sup{t > 0, doo(z n (t), dU) > CX}. 

We apply this lemma to the contagion model described in Section [A,3| Let us define, for 
< t < A 







_ y: s ^\t 

0<l<6 


«5-(r) 


j,k,e 


— t, and 


S(t) 





with s^' k ' e ' 1 given in Lemma A. 8. With Bin(j, ir) denoting a binomial variable with parameters 
j and 7r, we have 



5 j ' k ' 9 (T) 
<T(r) 



fx(j,k)p(j,k,9)¥[Bm(j,-)>e 

j,k,0 



T — T 



j,k,9<j 

T T 

X(I( T 



and 



f(r) 



i,fc,o<e<j 

A. 5 Proof of Theorem 13.61 



T T 

E Kj,k)p(j,k,e)¥(Bm(j,^)>e) . 



(16) 



(17) 



(18) 



We now proceed to the proof of Theorem 3.6 whose aim is to approximate the value D n (T n )/n 
as n — > oo. We base the proof on Theorem A. 9 However, several difficulties arise since in our 
case since the number of variables depends on n. We first need to bound the contribution of 
higher order terms in the infinite sums (|17|) and (|18[). Fix e > 0. By Condition 3.1 we know 



A = E Mi> k ) = E J/ ^' fc ) G (°> 00 ) 



Then, there exists an integer K e , such that 



E E^> fc ) + E ^2^U,k) <e, 

k>K E j j>K E k 
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which implies that 

jAk>K E 

It follows that 



V < r < A, fe A*0% k)pU, k, 9)¥ (Bin(i, I) > e) < e. (19) 



jAk>K s fl<6<j 



The number of vertices with degree (j, k) is n/j, n (j, k). Again, by Condition 3.1 



Y k Vn(j,k) = ^j/%0',fc) -)■ A G (0,oo). 
Therefore, for n large enough, Yl-jAk>K ^Mn(j>^) < £ 5 and 



V0<i<m n , ^ kD^ e (t)/n<e. (20) 

jAk>K e ,0<6<j 



For if > 1, we denote 



y* := (bSJ>M(t)) and 

V / jAk<K, 0<l<d<j 



jAk<K, 0<K6<j 

both of dimension b(K), where <5 3 ''*' (t) , s 3 "**' ' (t) are solutions to a system (DE) of ordinary 
differential equations. Let 

vr* = min{vr G [0, 1]|J(tt) = vr}. 
For the arbitrary constant e > we fixed above, we define the domain U e as 

C/ E = {( r , y K ') e R b ^ +1 : _ e < r < A - e , -e < s j ' k ' 6 ' 1 < 1}. (21) 

The domain U e is a bounded open set which contains the support of all initial values of the 
variables. Each variable is bounded by a constant times n (Co = 1). By the definition of our 
process, the Boundedness condition is satisfied with f3(n) = 1. The second condition of the 
theorem is satisfied by some Ai(re) = 0(l/n). Finally the Lipschitz property is also satisfied 
since A — t is bounded away from zero. Then by Lemma A. 9 and by using Lemma A. 6 for 
convergence of initial conditions, we have : 

Corollary A. 10. For a sufficiently large constant C , we have 

P(Vt < na c {n),Y% e {t) = ny Ke {t/n) + 0(n 3 / 4 )) = 1 - 0(6(K £ )n- 1 / 4 exp(-n~ 1 / 4 )) (22) 
uniformly for all t < no~c( n ) where 

a c {n) = sup{r > 0, d(y Ke (r), dU 6 ) > Cn" 1 / 4 }. 
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When the solution reaches the boundary of U e , it violates the first constraint in 21 determined 
by f = A — e. By convergence of — to A, there is a value no such that Vn > no, ^ > A — e, 
which ensures that fn < m n . Using (|19|) and (|20|), we have, for < t < nf and n > no: 



\D-(t)/n-S-(t/n)\ 



k(D^ e (t)/n-5^(t/n))\ 



kfii 



< 



E 



E 

jAk<K s 0<j 



k\D^ e (t)/n-5 j ' k ' 9 (t/n) 



+ 2e, 



and 



\D n (t)/n-5(t/n)\ < 



E 



E 



D J ' h 



\t)/n-5 j ' k ' e (t/n) 



+ 2e, 



We obtain by Corollary A. 10 that 



sup \D~(t)/n- 5~(t/n)\ < 2e + o p (l) 

t<fn 

sup \D n (t)/n - 5{t/n)\ < 2e + o p (l) 

t<fn 



(23) 

(24) 

(25) 
(26) 



We nw study the stopping time T n defined in (13) and the size of the default cascade D n (T n ). 
First assume I(tt) > ir for all ir £ [0, 1), i.e., ir* = 1. Then we have 



Vr < T,<r(r) = ^k5 3 



kfi, 



t) -t > 0. 



We have then that T n /n = f + 0(e) + o p (l) and from convergence (26), since 5(f) = 1 — 0(e), 
we obtain by tending e to that \D n (T n )\ = n — o p (n). This proves the first part of the theorem. 

Now consider the case tt* < 1, and furthermore tt* is a stable fixed point of I (it). Then 
by definition of tt* and by using the fact that 1(1) < 1, we have /(vr) < tt for some interval 
(tt*,tt* + tt). Then <E (r) is negative in an interval (t*,t* + f), with r* = Xtt* . 

Let e such that 2e < — inf rg ( r » iT *+f) $~ ( T ) an d denote a the first iteration at which it 
reaches the minimum. Since (E(<r) < —2e it follows that with high probability D~ (an) / n < 0, 
so T n /n = r* + 0(e) + o p (l). The conclusion follows by taking the limit e — > 0. 



A.6 Proof of Theorem Q 

Strong connectivity sparse random directed graphs with prescribed degree sequence has been 
studied by Cooper and Frieze in [18] . Let A n represent the average degree (then by Condi- 
3.1, A n — > A as n — > oo), and n n (j,k) represent the empirical distribution of the degrees, 



tion 



assumed to be proper (as defined below), then [181 Theorem 1.2] states that if 

> i, 

j,k 



(27) 
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then the graph contains w.h.p. a strongly connected giant component. 

We remark that the theorem above is given in [18J under stronger assumptions on the degree 
sequence, adding to Assumption 3.1 the following three conditions, in which A n denotes the 
maximum degree: 

• Let p n = max(£ M - £. ^ jH ^' j) ). If A n -> oo with n then p n = o(A n ). 

,1/12 



A„ < -j . 

n — logn 

As n — > oo, v n — ¥ y G (0, oo) 



Following |18j we call a degree sequence proper if it satisfies Assumption 3.1 together with the 
above conditions. 

A first reason for adding these conditions in [18] is to ensure that Equation (11) holds. 
However, following Janson [30], the restricted set of conditions 3.1 is sufficient. The second 
reason is that [18] gives more precise results on the structure of the giant component. For our 
purpose, to find the sufficient condition for the existence of strongly connected giant component, 
we show that these supplementary conditions may be dropped. 

It is easy to see that a bounded degree sequence (i.e., A n = O(l)) which satisfies Assumption 



3.1 is proper. We use this fact in the following. 



Lemma A. 11. Consider the random directed graph G*(d n ,d+) constructed by configuration 
model, where the degree sequence satisfies Assumption 3.1, If 

£;^>i, (28) 

then with high probability the graph contains a strongly connected giant component. 

Proof. By the second moment property and Fatou's lemma, there exists a constant C such that 

j,k j,k 

< liminfV(i 2 + A; 2 K(i,A;)<C. 
n— >oo L — * 

Then, it follows that for arbitrary e > 0, there exists a constant A £ such that 

jkfi(j,k)<e. 

jAk>A £ 

Thus, by choosing e small enough, there exists a constant A E such that 

jr'Afc<A e 

We now modify the graph such that the maximum degree is equal to A e : for every node i such 
that d+(i) A d~(i) > A e , all its in-coming (resp. out-going) half-edges are transferred to new 
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nodes with degree (0,1) (resp. with degree (1,0)). Since these newly created nodes cannot 
be part of any strongly connected component, it follows that, if the modified graph contains 
such a component, then necessarily the initial graph also does. It is then enough to evaluate 



Equation (27) for this modified graph, which by construction verifies the Assumption 3.1 for 
the new empirical distribution \x with the average degree A. Also, since the degrees of the 
modified graph are bounded, the supplementary conditions above also hold, i.e., the degree 



sequence is proper, and we can apply Cooper & Frieze's result. It only remains to show that 

A 



^/j k > 1. Indeed, we have 



j,k jAk<A e 



0<i,fc<A £ 



A 

p>(j, fc) 



0<j,k<A e 

The last equality follows from the fact that for < j, k < A e , we have 

AC/, fc ) = Mji fc) 
A A • 

This is true since the total number of edges, and the number of nodes with degree j, k for 
< j, k < A e , stays unmodified. 

□ 



We now proceed to the proof of Theorem 3.9 Our proof is based on ideas applied in [2TJ [29 



for site and bond percolation in configuration model. Our aim is to show that the skeleton of 
contagious links in the random financial network is still described by configuration model, with 



a degree sequence verifying Assumptions 3.1, and then apply Lemma A. 11 



For each node i, the set of contagious out-going edges is given by 

C n (i) :={l | (1-R)e n (i,l)> ln (i)}. 
Let us denote their number by 

4(*):=#C n (*). 

We denote by G c n the unweighted skeleton of contag ious links in the random network Cr,^(G n ), 
endowed with the capital ratios j n . 

In order to characterize the law of G^j, we adapt Janson's method [29J for the directed case. 

Lemma A. 12. The unweighted skeleton of contagious links has the same law as the random 
graph constructed as follows: 
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1. Replace the degree sequence (d+,d n ) of size n by the degree sequence (d+,,d n ,) of size n' , 
with 



n = n 



+ m n - <£(»)» 



i=i 



V 1 < i < n, d+,(i) = c+(i), d n ,{i) = d n (i), 
V n + 1 < i < n, d+(i) = 1, d+ (i) = 0. 

12. Construct the random unweighted graph G* n , (d^, , d~, ) wii/i n' nodes, and the degree se- 
quence (d^,,d~,) by configuration model. 

3. Delete n + = n' — n randomly chosen nodes with out-degree 1 and in-degree 0. 

Proof. The skeleton can be obtained in a two-step procedure. First, disconnect all non- 
contagious links in G*(e n ) from their end nodes and transfer them to newly created nodes 
of degree (1,0). Then delete all new nodes and their incident edges. The first step of this 
procedure may be dubbed as "rewiring". Looking at graphs as configurations, and since the 
first step changes the total number of nodes but not the number of half-edges, it is easy to 
see that there is a one to one correspondence between the configurations before and after the 
"rewiring". Thus, the graph after rewiring is still described by the configuration model, and 
has the same law as G* n , (d^", , d~, ) . Finally, by symmetry, the nodes with out-degree 1 and 
in-degree are equivalent, so one may remove randomly the appropriate number of them. □ 



Note that since the degree sequence before rewiring verifies Condition 3.1, so does the 
degree sequence after rewiring. Moreover, since we are interested in the strongly connected 
component and nodes of degrees (1,0) will not be included, we can actually apply Lemma 



A. 11 to the random graph resulting by the above contagion process. Hence, we may study the 
strongly connected component in the intermediate graph G*,(d^",, d~,). 

Let us denote by l n i(j,k), the number of nodes with out-degree j and in-degree k in the 
graph G* n , (d+, , d^, ) , and by \ n >, the average degree. Then the average directed degree in this 
random graph is given by v n := J2j,kj kl n'U^ k )/(^n'n'). 

We first observe that \ n m' = An, since the number of edges is unchanged after rewiring of 
the links. For every k > 0, the quantity ^ ■ jl n '(j, k) represents the number of out-going edges 
belonging to nodes with in-degree k in the graph after rewiring, which in turn represents the 
number of contagious out-going edges belonging to nodes with in-degree k in the graph before 
rewiring. But so does ^ - p n (j, k,l)nfj, n (j, k)j . So, for all k 

j j 

= 2>u,*,i)^ 

J 



A 



39 



where convergence holds by the second moment property in Assumption 3.1 Applying Lemma 

shows t 



A. 11 to the sequence of degrees in the graph after rewiring shows that when 

k j 



X n m' 



k j 



A 



> 1, 



then with high probability there exists a giant strongly connected component in the skeleton 
of contagious links. 
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